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We propose an approach to reconstruct any superconducting charge qubit state by using quantum state to- 
mography. This procedure requires a series of measurements on a large enough number of identically prepared 
copies of the quantum system. The experimental feasibility of this procedure is explained and the time scales 
for different quantum operations are estimated according to experimentally accessible parameters. Based on the 
state tomography, we also investigate the possibility of the process tomography. 
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I. INTRODUCTION 

The generation of superpositions of macroscopic quantum 
states in superconducting devices have moti- 

vated further research on quantum information processing in 
these systems. Two types of superconducting qubits based on 
Josephson junction devices have been proposed and experi- 
mentally demonstrated. One involves two Cooper-pair charge 
states in a small superconducting island connected to a cir- 
cuit by a Josephson tunnel junction and a gate capacitor (see, 
e.g., OSS])- An alternative approach is based on the phase 
states of a Josephson junction or the flux states in a ring su- 
perconducting structure 0.S01- Further, experimental ob- 
servations on quantum oscillations and the demonstration of 
conditional gate operations in two coupled charge qubits |3] 
are necessary first steps towards future realizations of quan- 
tum information processors. 

A crucial step in quantum information processing is the 
measurement of the output quantum states. However, a quan- 
tum state cannot be ascertained by a single quantum measure- 
ment. This is because quantum states may comprise many 
complementary features which cannot be measured simulta- 
neously and precisely due to uncertainty relations. However, 
all complementary aspects can in principle be observed by a 
series of measurements on a large enough number of identi- 
cally prepared copies of the quantum system. Then we can 
reconstruct a quantum state from such a complete set of mea- 
surements of system observables (i.e., the quorum |8]). Such 
a procedure is called "reconstruction of quantum states" or 
Quantum State Tomography (QST). 

Quantum state tomography is not only important for quan- 
tum computation, which requires the verification of the accu- 
racy of quantum operations, but it is also important for fun- 
damental physics. Many theoretical studies for tomographic 
reconstruction of quantum states have been done, e.g. refer- 
ences £3, Q3, [01 fl^l . Experimentally, tomography has been 
investigated for a variety of systems, including: e.g., the vi- 
brational state of molecules 1 13], the motional quantum state 
of a trapped atom fl4lfl5ll . two-photon states lift , the elec- 
tromagnetic field 11711 . and rare-earth-metal-ion-based solid- 



state qubit 1 18], the two-qubit states in the trapped ions fl9il . 
The quantum states of multiple spin-i nuclei have also been 
measured in the high-temperature regime using NMR tech- 
niques Enl EHEl . 

For continuous variable cases (e.g., the molecular vibra- 
tional mode fiUl . motional quantum states of a trapped 
ion 1 14, 11511 . a single-mode El of the electromagnetic field), 
the quantum states can be known by the tomographic mea- 
surement of their Wigner function. For the discrete variable 
case (e.g., in NMR systems), the measurements on the density 
matrix in NMR experiments are realized by the NMR spec- 
trum of the linear combinations of "product operators", i.e. 
products of the usual angular momentum operators l22Tl . 

Based on the state tomography, a quantum "black box" con- 
nected to an unknown external reservoir can also be charac- 
terized. This "black box" transfers any known input state to 
an unknown output state. The determination of the quantum 
transfer function for this 'black box" is called quantum pro- 
cess tomography 12311 . This procedure needs to input a large 
enough number of different known states into the "black box", 
then to make tomographic measurements on output states, fi- 
nally to obtain the quantum transfer function, which deter- 
mines the "black box". This procedure would be very impor- 
tant for the case when the noisy channel is unclear. Process 
tomography has been experimentally realized, e.g., in optical 
systems KM. NMR El. 

To our knowledge, there is no adequate theoretical anal- 
ysis or experimental demonstration for the reconstruction of 
qubit states in solid state systems, besides our recent work in 
Ref. 12611 . There, we considered a very general class of spin 
Hamiltonians used to model generic solid state systems 1 26] . 
Here, the emphasis is not on a general model but on a specific 
system: superconducting qubits. Recent technical progress 
makes it possible to realize quantum control in superconduct- 
ing quantum devices and ascertain either the charge | \ H or 
the flux |4] qubit states. Furthermore, practical experiments 
on quantum computing require the knowledge of the full in- 
formation of the quantum state, so the reconstruction of quan- 
tum states in solid state systems is a very important issue. 

In this paper, we analyze how to reconstruct charge qubit 
states in superconducting circuits. In principle, if all qubits 
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can be measured at same time in superconducting circuits as 
in optical systems (e.g., Ref. 1 10]), then only single-qubit op- 
erations are enough to assist the implementation of the re- 
construction of any multiple-qubit state. However, our pro- 
posal only considers one qubit measurement at a time. This is 
because simultaneous measurements of many qubits are cur- 
rently very difficult to implement l27il in superconducting cir- 
cuits. Another reason is that simultaneous measurements re- 
quire many probes in contact with qubits, inducing more noisy 
channels. These multiple noisy channels will quickly reduce 
the coherence of the qubit states, decreasing the accuracy of 
the reconstructed states. So for two-qubit and multi-qubit state 
tomography, appropriate two-qubit operations are necessary, 
due to the constraint of a single-qubit measurement at a time. 

Although our analysis of the tomographic reconstruction of 
charge qubit states might seem somewhat similar to the one 
used for NMR systems there are significant dif- 

ferences on how to realize the state tomography in Josephson 
junction (JJ) charge qubits. For example, NMR QST (like 
optical QST) also only involves single-qubit operations. A 
question we will focus on is the following: is it possible to 
do QST with the currently accessible experimental capability 
on J J qubits? In view of the short relaxation and decoherence 
times, it is also necessary to estimate quantum operation times 
required for reconstructing charge qubit states. In particular, 
it is not trivial to find an appropriate two-qubit operation to 
realize all two and multiple qubit measurements. 

Here, we theoretically analyze in detail the necessary ex- 
perimental steps for the tomographic reconstruction of dc- 
SQUID-controlled charge qubit states. This analysis can be 
easily generalized to other proposals of controllable super- 
conducting qubits (e.g., flux and phase qubits). In Sec. II, 
the reconstruction of single-qubit states is described in de- 
tail. The time scales of operations for measurements of all 
three unknown matrix elements are also estimated by using 
currently accessible experimental parameters. In Sec. Ill, all 
operations required to reconstruct two-qubit states are given, 
the time scales for the first and second qubit measurements 
are estimated using experimentally accessible parameters. In 
Sec. IV, using an example, we generalize our two-qubit to- 
mography to the multiple-qubit case. Finally in Sec. V, we 
discuss the "process tomography" of singe-qubit charge sys- 
tems based the "state tomography". Sections III, IV, and V 
contain our most important results. The conclusions and fur- 
ther discussions are given in sections VI and VII, respectively. 



II. RECONSTRUCTION OF SINGLE-QUBIT STATES 

The content of this section on single-qubit operations and 
the reconstruction of single-qubit states is known J2ql to spe- 
cialists in the optical, NMR and other areas (e.g, M@ ll(fllril 
d d 0O (S El E d IH IH H), where the QST is 
extensively studied. But here we specify a detailed descrip- 
tion of the steps needed for the experimental realization of 
the tomographic reconstruction for charge qubit states. This 
should be helpful to solid state experimentalists who are not 
specialists on the QST. 



A. Theoretical model and single-qubit states 

We consider a controllable dc-SQUID system which con- 
sists of a small superconducting island with n excess Cooper- 
pair charges, connected by two nominally-identical ultra- 
small Josephson junctions; each having capacitance Cj and 
Josephson coupling energy Ej. A control gate voltage V g is 
coupled to the Cooper-pair island by a gate capacitance C g . 
The qubit is assumed to work in the charge regime, e.g., the 
single-electron charging energy Eq = e 2 /2(C g + 2Cj) and 
Josephson coupling energy E® satisfy the condition Eq ^> 
Ej. If the applied gate voltage range V g is near a value 
V g = e/C g , only two charge states, denoted by n = and 
n = 1, play a key role, then this charged box is reduced to 
a two-level system (qubit) whose dynamical evolution is gov- 
erned by the Hamiltonian la. 12911 

H=-^5E ch (n e )a z -^Ej(<t> x )a x , (1) 

where we adopt the convention of charge states |0) = 1 1) and 
|1) = ||). The charge energy <5E ch (n g ) = AE G (\-2n g ) with 
n g = C g Vg/2e can be controlled by the gate voltage V g . The 
Josephson coupling energy Ej(Q x ) = 2 E® cos(ir$> x /&o) is 
adjustable by the external flux and <&o = h/2e is the flux 
quantum. Our goal here is to determine any single charge 
qubit state by the controllable dynamical operation governed 
by the Hamiltonian Q. 

Any single-qubit state (mixed or pure) can be represented 
by a density matrix operator in a basis { 1 0) = |T), |1) = \i)} 
as 




P = Poo|0)(0|+p O i|0)(l|+pio|l)(0| 

+ P11 11) <1|, (2b) 

where <Jk=x,y,z are Pauli operators and ak=o is an identity 
operator. Four real parameters r& (fc = 0, x, y, z) can be ex- 
pressed as 

ro = Poo + r x = p l + Pw, 

r v = i(p i - pw), r z = p 00 - p n . 

The normalization condition poo + pn = 1 ensures that the 
qubit (0 can actually be determined by three real parameters 
t x , t v , r z corresponding |28] to a Bloch vector r , which sat- 
isfies the condition | r | < 1 (see FigQa)). The state p is pure 
if and only if IT*! = 1. When the state p is pure, the Bloch 
vector T* defines a point on the unit three-dimensional sphere. 

These three coefficients rk (k — x, y, z) can be obtained 
from measurements of a x , a y , a z . The correspondence be- 
tween these three measurements and the coefficients r-^ is 
given by 

r k = Tr(p at), 

due to the relation Tr(cTjCrj) = 2<Jy, where 5^ is the Kro- 
necker delta. 



Quantum operations and measurements on single-qubit 
states 



In principle, the state of the charge qubit can be read by a 
single-electron transistor (SET) |2, 3, 29] coupled capacitively 
to a charge qubit. Here we consider the ideal case in which 
the SET is coupled to the qubit only during the measurement. 
When the SET is coupled to the qubit, the dissipative current 
/ flowing through the SET is proportional to the probability 
of a projective operator measurement 1 1) (1 1 on the qubit state, 
which has actually been applied by the experiment yjyj]. The 
|1)(1| measurement is equivalent to a a z measurement on the 
state p, 



Pi 



Tr(p|l){l|) 



~[1 -Tr(pcr 2 



Pn 



due to the relation 



|1>(1| = ^-<T Z ). 

The parameters tq and r z can be determined by the result of 
the measurement |1)(1|, together with the normalization con- 
dition. 

We can also relate the two other measurement operators, 
<j x and <j y , to the operator |1)(1| (essentially a z ), which is 
the measurement experimentally realized in the charge qubits. 
This is because the current / flowing through the SET is sen- 
sitive to the charge state |1), so the single qubit operations 
have to be performed so that the desired parameter r x or r y is 
transformed to the measured diagonal positions. 

Now, we describe the steps to measure r x or r y . Let us first 
choose the external flux <fr x = and suddenly drive the qubit 
to the degeneracy point for a time 



Htt 



2£j(0) 4E° 



such that the qubit state can be rotated — tt/2 along the x di- 
rection, here Ej(0) = Ej(® x = 0). 

The probability p2 of the measurement 1 1) (1 1 on the rotated 
state is 



P-2 



Tr (R x (t x ) pRt(t x )\l)(l\) 

Tr (exp ji^CTzj p expj-i^e^j |1)(1| 

Tr (p expj-i^j exp{«^(7 x } 



2^ + r y ), 



(3) 



where R x (t x ) = exp {iEj(0)a x t x /2h}, Eq. (|3} means that 
the measurement |1) (1| on the state rotated — tt/2 along the x 
direction is equivalent to the measurement a y , and the rotation 
—tt/2 of the qubit is equivalent to an inverse rotation of the 
measuring instrument, see Fig. [2 

In order to make the third measurement <r x , the qubit state 
needs now to be rotated —tt/2 (or tt/2) along the y direction. 
This can be done (e.g., — tt/2 rotation) as follows: 




FIG. 1 : (Color on line) The black Bloch vector indicates a qubit state; 
the (yellow arrows) r x , r y and r z represent the three components of 
the Bloch vector along the x, y and z axes. The and 1 in the north 
and south poles of the Bloch sphere denote the measured states |0) 
and |1), respectively. The measurement instrument is attached to a 
pole (e.g., "1") of the sphere. A —tt/2 rotation of the qubit state 
along the x direction is equivalent to a tt/2 rotation of the measuring 
instrument along the x direction. 



(i) Set $j; = <&o /2 and n g = 0; let the system evolve a time 

t z ,i = fnr/8Ec such that a rotation of — 7r/2 along the 
z direction is realized. 

(ii) After the time t Zt ±, set $ x = and n g = 1/2 and let 
the system evolve a time period t Xt i = Htt/2Ej(0) = 
Htt/AEj such that the system rotates — 7r/2 along the x 
direction. 

(iii) Set & x = and n g = 1/2 again and let the system 
evolve a time t z ^ — 3IItt/8Ec and a rotation — 37r/2 
along the z direction is obtained. 

Combining the above three steps, shown in Fig. [2] a —tt/2 
rotation of the qubit along the y direction is realized. 

(iv) After the above rotations, a measurement 1 1) (1 1 on this 
rotated state must be made, which is equivalent to mea- 
suring a x . Then, the measured probability becomes 

P3 = Tr(R w pRl }XjZ 

= ^(exp{i~^}pexp{-i~^,}|l)(l|J 

= \0- + r*) 



RSz,l)Rx(t x ,l)R z (t z ,2),Wd 

2E C 



with R z . x . z 

Rz(tz,l) = 

Rx(t x ,i) = exp |i^-cr x t X) i| = exp|i^(T x | 



exp ^ t—^-a z t zA 



exp 



Rz(t„, 2 ) 



t 2E C . 

exp < l - (7 z t z 2 



3tt 

exp <; i— a z 
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FIG. 3: Graphical representation of the density matrix p for single- 
qubit states, see the example explained in section II. The real p 



(R) 



and imaginary p\- parts of the density matrix elements pij — (i\p\j) 
(i, j = 0, 1) are plotted in (a) and (b), respectively. 



FIG. 2: (Color on line) The Bloch vector is the same as Fig. Q A 
— 7r/2 rotation of the qubit along the y direction is equivalently re- 
alized by the rotation n/2 of the measuring instrument along the z 
direction (from (a) to (b)), then a 7r/2 rotation along the x direction 
(from (b) to (c)), and a 37r/2 rotation along z direction (from (c) to 
(d)). 



We explained how to measure the single qubit states by sin- 
gle qubit operations and measuring |1)(1|. Below, we give 
an example that shows a reconstructed single-qubit state can 
be graphically represented, and we further give estimates of 
the operation times to obtain each of the matrix elements of 
single-qubit states. 



C. An example 

The three measurement results (pi , p2, P3) can be used to 
obtain four coefficients (ro, r x , r y , r z ) that define a single- 
qubit state. A single-qubit state can be reconstructed follow- 
ing the steps presented above and an example is described 
here. If we obtain r x = 1, r y = \/3, r z = 1 by the three 
experimentally measured probabilities (pi, pi and ^3) on a 
quantum ensemble of an unknown charge qubit state p, then 



Poo 



1 

Pxx = 2> 



An 



1 = \ (l-n/3 



1 



Pw = 2 • 1 + 



Thus, the reconstructed state p can be written as 
1 



P 



2 

1 r 

4 



(|0}(0| + |1)(1|)4 

(l-iVa) |o><i] 



l+iVa) ii)(oi 



whose real p^P and imaginary pr) parts are graphically rep- 
resented in Fig. [3] 



D. Operation time estimates 

The coherent operations required for the tomographic mea- 
surements are limited by the decoherence time T2- Now let 
us explore whether the single-qubit state can be reconstructed 
with the current experiments. To estimate the corresponding 
time scales for quantum operations to obtain the measure- 
ments of a„ and a x , we first take the suggested parameters 
from Ref. j^, that is, E$ = 100 mK (about 8.6 /xeV or 2.08 
GHz) and E c = 1 K (about 86 peW or 20.8 GHz). Here, we 
use temperature units for energies as in reference 1 29] . Thus 
the approximate time scales of one-qubit operations to obtain 
r v and r x are 



t x w 5.9 x 10~ n s 



and 



in 



t z ,i + t x .i + t z 



7.1 x 10- 



These time scales, required to reconstruct the single-qubit 
states, are within the measured values |2, 3] of the decoher- 
ence time T2 (of the order of magnitude of ns) of single-qubit 
charge states. 

Now let us consider another set of experimental values. For 
example, if the Josephson and charge energies are taken (sec- 
ond paper in Ref. y|) as 2E^ = 45/xeV (about 520 mK or 
10.9 GHz) and 4E C = 580peV (about 6.73 K or 140 GHz), 
then the time scales required to reconstruct single-qubit states 



are about t 7 



2.3 x 10 _il s and t 



v 



3.0 x 10 _li s, which 



are within the decoherence time T2 = 5 ns obtained by that 
experiment y[ . 

If we take the Josephson and charge energies from in 
Ref. 03], that is, 2E°/h = 13.0 GHz (about 625 mK or 53.7 
/zeV) and 4E c /h = 149.1 GHz (about 7.16 K or 618 peV), 
then the time scales required to reconstruct single-qubit states 
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are about t x sa 1.9 x 10~ n s and t y w 2.6 x 10~ n s, which 
are also less than one order of magnitude of the decoherence 
time T 2 = 325 ps measured by that experiment [30]. 



and = (2C°)~ 1 + C^ 1 . Thus, the interaction between 
the two qubits can be controlled by two external fluxes <&i x 
applied to each qubit. 



III. RECONSTRUCTION OF TWO-QUBIT STATES 



B. Quantum operations and measurements on two-qubit states 



A. Theoretical model and two-qubit states 

In this section, we focus on the reconstruction of two-qubit 
charge states. Any two-qubit state pi can be characterized by 
a density matrix operator 



Pi 



(J\i <K> cr 2 j 



(4) 



i,j—0,x,y,z 



where the 16 parameters r£j are real numbers. The normal- 
ization property of the quantum state requires that ro,o = 1> 
so the state p\ in Eq. (0} can in principle be reconstructed 1I3U1 
by 15 measurements described by the operators an £g> a 2 j, 
where all i and j are not simultaneously taken to be 0. If one 
of <7ij (i = 0, x, y, z) or a 2 j (i = 0, x, y, z) is an iden- 
tity operator among the measurement operators a\i ® a 2 j, we 
call such a measurement a single-qubit measurement and only 
write out the non-identity Pauli operator in the following ex- 
pressions. For example, the operator a\ x ® 020 is called a x 
measurement of the first qubit, and abbreviated by <j\ x . So 
there are nine two-qubit measurements among measurements 
&U <8> <J2j • Recall that only one-qubit is involved during the 
measurement process in our approach. If we want to obtain 
these nine two-qubit measurements, then two-qubit operations 
must be applied |32] such that the singe-qubit measurement 
can be equivalently transformed into expected two-qubit mea- 
surements. 

Now our task is to find a non-local two-qubit operation and 
use this operation to realize all necessary two-qubit measure- 
ments on two-qubit states. Here we consider a model pro- 
posed by Makhlin et al. ||6|], where two charge qubits are cou- 
pled in parallel to a common inductor with inductance L. The 
Hamiltonian |0] is 

1 2 

H = --22[SE ch (ni ig )ai z + Ej(<S> lx )(T lx ] 



1=1 



E in t($lx,$2x)viy <3<J 2y , 



(5) 



where it is assumed that both qubits are identical, so the 
charge energies SE^n^g) and Josephson coupling energies 
Ej($i x ) take the same form as in Eq. 0, but now 8E c h(ni ig ) 
and Ej($>i x ) for each qubit can be separately controlled by 
the gate voltages and external fluxes. The interaction energy 
Ei nt for two coupled qubits can be written as 



E in t(<f>lx,®2x) 



E\, 



with 



E v 



Now, we discuss how to reconstruct two-qubit states from 
the experimental measurements (|1)(1|)/ (I = 1, 2). Single 
charge qubit operations can be realized by controlling the gate 
voltage and Josephson couplings. However the two-qubit op- 
erations need to couple a pair of interacting charge qubits. The 
realization of the coupling of two charge qubits have to simul- 
taneously turn on the Josephson couplings of the two charge 
qubits in Eq. then ai x terms have to be included in the 
two-qubit operation. However the charge energies for two 
qubits can be switched off by applying gate voltages such that 
7i; ig = 1/2 (I = 1, 2), so a two-qubit operation can be gov- 
erned by a simpler Hamiltonian 

H ' = -\Yl vix-E in t(<S> lx ,<S> 2x ) a ly <E>(T 2y , (6) 

1=1,2 

where charging energies are set to zero, SE^n^g) — (/ = 
1, 2), with ni g = 7i2 >g = 1/2, and the external magnetic 
fields are chosen such that 

*ix = $2* + |(2g+l)$ 

with the positive integer number q. The coupling 
E- m t($ix, $2x) can be controlled by the external fluxes $1^ 
and <& 2x . 

The basic two-qubit operation can be given by the time- 
evolution operator U(t) = cxp{—iH't/h}, which can be 
written by using the Pauli operators as 

1 sin 6' 

u (t) = 2 ( cos ^' + cos6') I + in z —^- {a lx + a 2x ) 

sin <// — n x sin 9' 

+ i ^ 0\z ® ®2z 



+ 1- 



sm + n x sm ( 



■ <J\y <X> cr 2y 



cos — cos ( 



' 0"1 X <£> <7 2x , 



(7) 



where 



1 



Ej 



?' = -£ int ($ lx ,$2x) Vl + a 2 . 
n 



Since the two external fluxes satisfy the condition <&i x = &2x, 
we let Ej{<5>\ x ) = Ej(§ 2x ) — Ej in the expression Eq. @ 
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for the two-qubit operation. The physical meaning of the an- 
gles 8' and <j>' becomes clearer by virtue of the "conjugation- 
by-^S" operation l33ll on the time evolution operator U(t), 
which is defined as 

U'(t) = exp |^((7i 1/ + a 2y )j U{t)exp{-i^(a ly + a 2y )} , 

here, £ = o\ y + a 2y . In the conjugate representation, the 
time evolution U'(t) corresponds to rotations ll34ll around the 
y axis by an angle <fi' and the (n x , 0, n z ) axis by an angle 8'. 
By choosing the duration t and tuning the values of Ej and 
Eint(&ix, $2x), we can obtain any desired two-qubit opera- 
tion. 

From Eq. 0, it is known that six single-qubit measure- 
ments {an, <J2j} with i, j — x,y,z and nine two-qubit mea- 
surements {an <E) <T2j} with i,j — x, y, z are enough to 
obtain fifteen parameters of the two-qubit state p\ = 
\ Hi,j=o,x, y ,z r hi ® a 2j- The single-qubit measurements 
(|l)(l|)i = \{<Jin - <Jiz) (I = 1, 2) on a given state pi can be 
obtained as follows. Two single-qubit measurements ai z and 
a 2z can be implemented by the direct measurements (|l)(l|)i 
and (|1)(1|)2 on the given state pi. Other four single-qubit 
measurements (corresponding to a% x , a 2x , &iy, (T 2y ) need 
single-qubit operations. 

The single-qubit operations corresponding to measure- 
ments <T[ X and <j\y on two-qubit states are the same as measur- 
ing a x and a y on single-qubit states. However, in the single 
qubit operations, we need to switch off the interaction of the 
two qubits. For example, in order to obtain the measurement 
<j\y, we need to switch off the interaction between the two- 
qubit system by setting the applied external flux $22: = tc/2, 
and setting the first subsystem at the degeneracy point and 
evolving a time t = Htt/4Ej. Finally, we make a measure- 
ment (|1)(1|) 1 on the rotated state, then the coefficient r y 
can be obtained by this measured result. The other three mea- 
surements can also be obtained by taking single-qubit opera- 
tions similar to cr ly . 

The single-qubit measurements have been obtained by the 
measurements (|1) — 1,2) on given states by using ap- 

propriate single-qubit operations as described above. In or- 
der to find out how to obtain the two-qubit measurements via 
(|l)(l|)z, let us consider the measurements (|1)(1|)/ on the 
given state p\ performed by a sequence W of single-qubit and 
two-qubit operations. The corresponding measured probabil- 
ity p can be expressed as 



P 



Tr[W Pl WtfllXll), 



-Tr [ Pl W^a lz W] , (8) 



where we show that the measurement (|1)(1|); on the rotated 
state W pi W' may be interpreted as an equivalent measure- 
ment W^aizW (I = 1, 2) on the state p\. So our task now 
is to find an appropriate two-qubit operation and apply this 
two-qubit and single-qubit operations to the measured state 
pi, such that we can equivalently obtain the desired two-qubit 
measurement. 

Here, the required two-qubit operation U(t) can be ob- 
tained by choosing the evolution time t, the Josephson cou- 
pling energies Ej, and E^ in Eq. Q such that <f>' = (2m — 



l)7r/4 and 8' = nn where to, n are positive integers. The 
above conditions can be satisfied if the ratio 




and the evolution time r is chosen as 

hn 



4£j 



V(4n) 2 - (2m -l) 2 



If we choose the integers m and n to minimize the ratio 
El/Ej, then E L /Ej = VTE = 3.87 when 8' = n and 
4>' = 7r/4; so the two-qubit operation time r is chosen as 
t = tm^/Tb/AEj. Thus Eq. is specified by the time evo- 
lution operator 



U(r) = 



1 



2V2 



1 - V2J I - (l + V2) (Tlx ® <T2x 
T2y +i<Jlz® C 2 z] • (9) 



Combined with other single qubit rotations, U(t) can be used 
to obtain all the desired coefficients r^j corresponding to the 
two-qubit measurements an <8 a 2 j, with i, j = x, y, z. 
Let us further discuss how to obtain a desired coefficient, 



for example, 



corresponding to the two-qubit measure- 



ment ai v ® a 2y . We can take the following steps: 

(i) We switch off the interaction between the first and sec- 
ond qubits by applying an external flux $2x = tt/2, 
which means £ , j($2x) = 0. Now we only manipu- 
late the first qubit such that a rotation tt/2 about the z 
axis, defined as Zi = exp[i7r<7i z /4], is performed; this 
single-qubit operation is described in Section II. 

(ii) Following the single-qubit rotation Z\ of the first qubit, 

the gate voltages are applied such that ni )g = — 
1/2, which means that the two qubits work at the de- 
generacy points. Simultaneously, we turn on and ad- 
just the external fluxes so that the external fluxes §i x , 
energies E^ and Ej in the two-qubit operation de- 
scribed by the Hamiltonian satisfy the conditions 
&i x = $2x 7^ 7r(2(j + l)/2 with positive integer q and 
El/Ej = y/TE = 3.87. Afterwards, we let the system 
evolve a time r = hn\/T5/4Ej; which means that a 
two-qubit rotation U (r) has been performed. 



The operation sequence W 
changes state pi into 



U(t)Zi described above 



p=U(r) Z x pi z\XJ\t). 

(iii) Finally, when a single-qubit measurement (|1)(1|) 1 is 
performed on the state p, a two-qubit measurement 
equivalent to ai z (8 a 2y is implemented: 



Z\ U\t) (|l)(l|) a U(t) Zi = - + — ^ (a u + a Xy ® a 2y ) 



The corresponding measurement probability p can be given as 
p = Tx{u{r)Z 1 p 1 Z\U\r){\l){l\) 1 ) 



1 



Tr[pi (a lz + a ly (g> tT 2y )} 



2V2 
1 , 

+ 2Vl + ' 



(10) 



Because the coefficient r 2j o = Tr(pi fiz), corresponding to 
the operator a\ z ® <720> nas been given by the single-qubit 
measurement a\ z , then the coefficient r y>y — Tr(pi u\ y ® 
<?2y) is obtained via p and r z .o- 

In table [I] we have summarized nine equivalent two-qubit 
measurements described by — \[2 W^a\ z W on the original 
state px, which are obtained by the first qubit measurement 
(|l)(l|)i on the rotated state W p\W^ for a sequence W of op- 
erations with appropriately-chosen single-qubit and two-qubit 
operations. We can use the results corresponding to these nine 
equivalent two-qubit measurements together with the other six 
single-qubit measurements to obtain all the coefficients corre- 
sponding to the two-qubit states, and then obtain any two qubit 
state. 

We can also obtain coefficients {i, j ^ 0) correspond- 
ing to all two-qubit measurements by using the second qubit 
measurement (|1)(1|)2. For example, if we make a measure- 
ment (| 1) (1|)2 on the rotated state p considered above, we 
obtain another equivalent two-qubit measurement, which is 
expressed as 



Z 1 t [/t(r)(|l)(l|) 2 C/(r)Z 1 



1 1 / 



Using this measurement, combined with the single-qubit mea- 
surement (T2 2 , we can obtain the coefficient r x<x correspond- 
ing to the two-qubit measurement <j\ x ® a^ x - Nine equiv- 
alent two-qubit measurements realized by the second qubit 
(|1)(1|)2 have also been summarized in Table HT1 Compar- 
ing tables U and |ll| shows that different operations and steps 
are required in order to obtain the same coefficient for dif- 
ferent measurements. For example, in order to obtain r x z , 
two operation steps are needed for the first qubit measurement 
( 1 1 ) ( 1 1 ) 1 , but it needs four steps for the second qubit measure- 
ment(|l)(l|) 2 . 



C. An example 

We can also give another schematic example for a recon- 
structed two-qubit state. For instance, according to the opera- 
tions steps discussed above for the reconstruction of any two- 
qubit state, if we obtain r XtX — 1/8, r XjV = r VtX — V3/8, 
t z ,z = l/4andr a!a = —1/8 from the sixteen measured prob- 
abilities on an ensemble of identically prepared copies of a 
two-qubit system with unknown state p', then we can recon- 
struct this unknown state as 



P 



|00)(0Q| + |ll)(U|] + i(l-iV3)|00)(ll| 
l + iVs) 1 11> <00|, 



p (R) 



(a) 



ij.kl 



(b) 





FIG. 4: Graphical representation of the density matrix p' for the two- 
qubit state described in the example given in section III. The real 
Pijkt ali d imaginary p'^ kl parts of the density matrix elements for 
the two-qubit state p' in the basis |00), |01), 1 10) , |11) are plotted 
in (a) and (b) respectively. 



which is graphically shown in Fig. @ with the real ~ 1 



and imaginary pfj kt parts of the reconstructed state p', where 
i, j, k, I can take the values or 1. 



Pij,kl 



D. Operation time estimates 

We can also estimate the operation time required to re- 
construct two-qubit states for the Josephson and charge en- 
ergies E° } = 100 mK and E c = 1 K. We assume that 
the ratio E^/Ej = \/l5 = 3.87 is obtained by adjusting 
the external flux &i x (I = 1,2) such that $i x = 0, which 
means the ratio between E^ and Ej should satisfy the con- 



dition E h /E^ = 2\/l5 = 7.74 when the circuits are fabri- 
cated. In such case, the realization of the two-qubit operation 
in Eq. l|9} requires a time r « 2.32 x 10~ 10 s. Our previous 
estimates for the times to perform ir/2 rotations about the x 
and z axes are 5.9 x 10~ n s and 3.0 x 10~ 12 s, respectively. 
Then, using tables|I|and|nj we can estimate the total operation 
time required for obtaining the coefficients of the two-qubit 
measurements corresponding to the first or second qubit mea- 
surements, respectively. We find that the required operation 
times for the two-qubit measurements are less than 0.4 ns for 
the two-qubit measurements. The decoherence time T2 (e.g., 
the decoherence time of charge qubit is about 5 ns in refer- 
ence 1 2]) experimentally obtained shows that it is possible to 
reconstruct two-qubit states within the current measurement 
technology. 

At present, completely controllable multi-qubit supercon- 
ducting circuits are not experimentally achievable. Here, let 
us consider the operation time estimates based on another con- 
trollable model 1 35]. In this model, N charge qubits are cou- 
pled to a common superconducting inductance L. The Hamil- 
tonian of any pair of qubits, say i and j, is 



= E ( £ *4 fc) + <*<4*°) + x4 l) ® 4°, (id 

k=i,j 



where the coupling constant x can be tuned to zero by chang- 
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TABLE I: Equivalent two-qubit measurements — y/2 o\ z W ob- 
tained by measuring (|l)(l|)i on the state W pi with a sequence 
of appropriately-chosen quantum operations W . 



TABLE II: Equivalent two-qubit measurements — \f2 W* 02 Z W ob- 
tained by measuring (j 1) { 1 1 ) 2 on the state W pi W' with a sequence 
of appropriately-chosen quantum operations W . 



Two-qubit 


Quantum 


Equivalent two-qubit 


Two-qubit 


Quantum 


Equivalent quantum 


measurement 


operation" W 


measurement 


measurement 


operation W 


measurement 


Tlx 6. 


5 02y 


U(t) 


Olz + Olx € 


) 02y 


Olx i 


5 02x 


U(t)Zi 


02z — Olx ® 02x 


Olx i 


5 02z 


XiU(t) 


— Oly + Olx ' 


S> 02z 


Oly e. 


5 2x 


U(r) 


02z + Oly C2> 02x 


Olx 6. 


5 02x 


U(t)Z 2 


Olz — Olx € 


) <X2ic 


Olz <5 


i 02x 


U(t)Xi 


— 02y + Olz ® 02x 


Uly 6. 


5 <T2y 


U(t)Zi 


Olz + Oly (§ 


)02y 


oix 6. 


5 02y 


U{t)ZiZ2 


02z — Olx <g> 02y 


<Jly 6. 


5 02z 


XiU(t)Zi 


Olx + Oly (§ 


) CT 2z 


Oly i 


<5 02y 


U(t)Z 2 


02z + Oly <g> 02y 


Oly £ 


5 02x 


U(t)ZiZ 2 


Olz — Oly g 


) O^z 


oiz 6. 


3 02y 


U(t)XiZ 2 


02x + Olz CS> 02y 


Viz S 


5 02y 


U(r)ZiXi 


— <Tl a + Olz ' 




Olx t 


5 0"2z 


U(t)ZiZ 2 X 2 


— 02y — Olx <S> 02z 


<Tlz £ 


5 <T2z 


XiU{t)ZiXi 


Olx + Olz € 


) 0" 2z 


Oly C. 


5 cr 2z 


U(t)Z 2 X 2 


— 02y + Oly ® 02z 


Viz <5 


5 02x 


U(t)ZiZ 2 Xi 


— Oly — Olz < 


S> 02x 


Olz i 


5 cr 2z 


U(r)XiZ 2 X 2 


02x + Olz ® 02z 



"Xi and Zi denote single qubit rotations 7r/2 of 2th qubit about the x and z 
axes, respectively, and r = Htt\/T5/4:Ej. 



ing the flux either through the common inductance L, or 
through the qubit i (or j). Moreover, the parameters Ek and 
£/t are respectively controlled by the voltage applied to the fcth 
qubit and the magnetic flux through the fcth qubit. The condi- 
tional logic gates, e.g., controlled-NOT and controlled-phase- 
shift gates, can be performed by virtue of only one two-bit 
operation and also single-qubit operations in this circuit. This 
approach is more accessible to experiments, facilitating to- 
mographic measurements. According to calculations 1 26] of 
tomographic measurements for a class of representative quan- 
tum computing models of solid state systems, the two-qubit 
operation required for the realization of the multi-qubit mea- 
surements in this circuit can be easily obtained. That is, if the 
ratio between the Josephson energy E® and the two-qubit cou- 
pling energy x is E® = 2x, when the circuit is fabricated, then 
a two-qubit operation U (t') = — io\ x ® <J2 X can be obtained 
with the evolution time r' ~ 1.2 x 1CP 10 s when the Joseph- 
son energy is taken as = 100 mK. Here, we assume that 
the two charge qubits are identical and the Josephson energies 
are maximum when the two-qubit operation is performed If 
the charging energy is taken as Eq = 1 K, then ir/2 rota- 
tions around the z and x axes need times 3.0 x 10~ 12 s and 
5.9 x 10~ n s, respectively. The operations to get each of the 
sixteen (single- and two-qubit) measurements can also be ob- 
tained for this model by using an approach similar to the one 
described above, the estimated operation times to obtain all 
coefficients of the two-qubit measurements are less than 0.3 
ns, which is also within the experimentally obtained decoher- 
ence time T2 = 5 ns. 



IV. RECONSTRUCTION OF MULTIPLE QUBIT STATES 

In the above two sections, we focused on the reconstruction 
of the single and two qubits states. In this section, we dis- 
cuss the reconstruction of any n-qubit state. In the multiple 
qubit charge circuit, the dynamical evolution is governed by 



the Hamiltonian I 



1 " 

H = - - ^[SE ch (n l:g ) a lz + Ej(<f> lx ) a lx ] 
2=1 

- ^ Ej nt ($l x ,§kx) Vly ®(Jky, 



(12) 



Kk 



where 5E ch {n Ug ) = AE C {1 - 2n; jg ), Ej($ lx ) = 
2E%cos(tt§i x /$o), and E int ($i x , § k x) take the same form 
as in Eq. 0. We also assume E L /2E^ = VIE = 3.87 and the 
single-qubits are nominally identical. By virtue of the control- 
lable Hamiltonian d!2i . in principle we can use (n — 1) two- 
qubit operations together with some single-qubit operations to 
reconstruct any n-qubit state, which can also be described by 
the density matrix operator 



Pi 



1 

2" 



1 1 ,^2 ,■ ' ' ,ln— Q,x,y,z 



with 2" real parameters ri y i 2 ... ,/„ corresponding to the mea- 
surements ai t (g) cri 2 ■ ■ ■ ® er;„ . But, here, we only show how to 
obtain a coefficient corresponding to a three-qubit measure- 
ment. The generalization to obtain coefficients of multiple 
qubit measurements is straightforward. 

In order to determine a three-qubit state, we need to make, 
single-qubit, two-qubit, and three-qubit measurements. It is 
known that all coefficients corresponding to single-qubit and 
two-qubit measurements can be obtained by using the same 
operations and measurements (|l)(l|)i=i, 23 as in section I 
and II. When we make two-qubit operations on, for exam- 
ple, the first and second qubits, the interaction of the third 
qubit with these two qubits is switched off by the applied flux 
$3x = ir/2. Now let us show how to obtain the coefficients 
corresponding to the three-qubit measurements. For example, 
for the coefficient r x ^ z ^ y of the measurement o\ x ® 022 <8> 03y, 
we should make the following sequence of quantum opera- 
tions: 

(i) Switch off the interaction of the third qubit with the first 
and second qubits by applying the flux $3^ = 7r/2. 
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Then make a two-qubit operation U%2 (r), with the same 
form as Eq. We use the subscript "12" to denote 
two-qubit operations on the first and second qubits. 

(ii) Switch off the interaction between the first and second 
qubits by setting Q 2x = 7r/2, and making a tt/2 rotation 
about the z axis for the first qubit. 

(iii) Make another two-qubit rotation C/ 13 (t) on the first and 
third qubits by adjusting the external fluxes such that 
$i x = $3^ = 0. The two-qubit operation U\ 3 (t) takes 
the same form as Eq. but the subscript "2" of the 
Pauli operators in Eq. l|9} is replaced by the subscript 
"3". This process can be described as 

, 2 Ul2 (r ) p2 ( T ) ^ Zl U 12 (r) p 2 U\ 2 (r) Z\ 

^ U 13 Z x U 12 (r) P2 U\ 2 (r) Z\ u\ z . (13) 

(iv) Finally, make a measurement ( 1 1) (1 1 ) i on the above ro- 
tated state, and obtain the equivalent measurement 

uUt) Z\ u\ 3 (|l)<l|)i U 13 Z x Uxi=\~ \°u + 
+ \ {crix <& <?2y + &ly <8 o 3y - a lx (g> o 2z (g) a 3y ), (14) 



(ii) After a certain time, the output states evolve into un- 
known states. 

(iii) By using the state tomography, we can ascertain these 
unknown states. 

(vi) Finally, an unknown quantum channel is determined by 
the data obtained for the tomographic measurements on 
these states. 

Experimentally, in order to determine the noisy channel of the 
studied iV-qubits l23ll . TV 2 known states need to be prepared, 
and these states must have density matrices which span the 
space of any allowed input state density matrices. 

We have shown that single-qubit state tomography is ex- 
perimentally accessible. In order to perform quantum process 
tomography for a single charge qubit. Four kinds of different 
charge states |0), |1), (|0) + |1))/V2, and (|0) + i|l))/V2 
need to be experimentally prepared. These states can be gen- 
erated in a SQUID-based charge qubit with current experi- 
ments I2l l3l 13011 . Thus, the process tomography of a single 
charge qubit is achievable using current technology. With fur- 
ther developments of this technique, the process tomography 
of multiple charge qubits could also be realized, when data 
from multi-qubit state tomography is obtainable. 



and corresponding measurement result p" is 

// _ 1 _ r z,0,0 + r x,y,0 + r y,0,y ~ r x,z,y 

P - J 4 ■ W 

Finally, we can obtain the coefficient r x ^ y ^ z based on p" and 
the single and two qubit measurement results r z ,o,o, ^a;,j/,o 
and r y ,o y y, which can be obtained by using the same way de- 
scribed in sections II and III. Other coefficients correspond- 
ing to three-qubit measurements can also be obtained by us- 
ing a similar procedure. According to the estimated time 
for reconstructing the two-qubit states, we believe that it is 
also possible to reconstruct the three-qubit states using cur- 
rent technology. Any multiple-qubit can also be reconstructed 
by sequentially designing the single-qubit and two-qubit op- 
erations. The generalization to multiple-qubit is an extension 
of the procedure that we outlined above. 



V. QUANTUM PROCESS TOMOGRAPHY 

It is worth briefly reviewing that, based on qubit state to- 
mography, the noisy channel (usually denoted as the "black 
box") of the controllable charge qubits can also be deter- 
mined. This experimental determination of the dynamics of 
the "black box" is called quantum process tomography 1 23], 
which can be described as follows: 

(i) Many known quantum states of the system under in- 
vestigation are input into the "black box", which is an 
unknown quantum channel, for example, an arbitrary 
environment. 



VI. CONCLUSIONS 

In conclusion, we discuss how to reconstruct charge qubit 
states via one-qubit measurements using controllable super- 
conducting quantum devices. Detailed operations for recon- 
structing single- and two-qubit states are presented. Any n- 
qubit state can also be reconstructed by using n — 1 two-qubit 
operations similar to Eq. ij^Ji for different qubit pairs and com- 
bining these with required single-qubit operations. Thus the 
non-local two-qubit operation Eq. (|9jl plays a key role in the 
reconstruction of the multiple-qubit states. However, this two- 
qubit operation is not unique for achieving our purpose. We 
should note that operations to obtain a fixed coefficient corre- 
sponding to multiple-qubit measurements are not unique. The 
measurements (|1)(1|)/ (I = 1, 2, . . . , ri) on the given state 
with fixed operations W are different for each qubit I, because 
W is not symmetric when exchanging /. Our proposal can also 
be generalized to other superconducting charge qubit circuits 
with the coupling mediated by photons or a tunable oscillator, 
e.g., Refs. ||3g,l33, or other types of superconducting qubits, 
e.g., Refs. SEallllilH. 

We find that the longest operation times to obtain the coef- 
ficient of single-qubit and two-qubit states are of the order of 
0.01 ns and 0.4 ns, respectively, which is less than the deco- 
herence time |2] T 2 = 5 ns. Moreover, the tt and tt/2 pulses 
for single-qubit operations can be performed very well, e.g., in 
the experiments of the charge echo 1 43 ] , and NMR-like exper- 
iments 1 44]. Another experimental estimate shows us that the 
manipulation accuracy can reach 80 — 90% (e.g., as in the sec- 
ond reference of Ref. [2]). Thus the single-qubit states could 
be reconstructed and the process tomography should also be 
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accessible in single-qubit charge systems with current experi- 
mental capabilities. In principle, the two-qubit states can also 
be reconstructed by virtue of well-controlled time for the two- 
qubit operation. We should also note that larger values of the 
charge energy E^, the Josephson energy E®, and coupling 
energy E- m t($>ix, &kx) can make the operation times shorter. 
Thus these larger values should be realized in order to facili- 
tate the tomographic reconstruction. 

Quantum oscillations and conditional gate operations have 
been demonstrated in two coupled charge qubits with the in- 
teractions 1 3] always turned on. Completely controllable two- 
qubit charge systems have not been realized yet. However, the 
coupled two charge qubits, allowing on and off switching of 
the interaction, might be realizable in the future ll45ll . Then 
our proposal will become realizable. Because the unswitch- 
able two-qubit interaction makes single-qubit operations im- 
possible, our proposed scheme cannot be readily used to the 
experimental reconstruction of multiple-qubit charge states 
when the two-qubit interactions are always turned on. How- 
ever, for the two-qubit circuit with "always-on" interaction, 
most of the single-qubit parameters l38t l39ll40H4lll42il can 
be tuned. We can adjust these parameters to obtain 15 dif- 
ferent two-qubit operations, and then derive 15 different mea- 
surement equations with these operations on input states. Af- 
terwards, the two-qubit states can finally be determined. The 
details on how to reconstruct the superconducting two-qubit 
states with the "always-on" couplings will be presented else- 
where. However, how to reconstruct qubit states in multiple- 
qubit (more than two qubits) circuits with "always-on" inter- 
actions is an open problem. 

VII. DISCUSSIONS 

In our paper, to simplify the algebra, we focus on one par- 
ticular family of measurements which are constructed by the 
direct product of the Pauli operators. However, one can con- 
ceive that other complete sets of measurements can also be 
used to do tomography. These different complete sets of mea- 
surements can be transformed to each other by unitary opera- 
tors. In practice, within the duration of the controllable manip- 
ulation, the smaller Bloch rotation might be advantageous to 
speed up the measurements, but it might also decrease the ac- 
curacies of the measurements due to a longer measuring time. 
How to choose suitable sets of operations during the measure- 
ment process is an important technical question for the mea- 
surement. 

It should also be pointed out that here we discuss an ideal 
case. In practice, the environmental effect is unavoidable, 
which result in the relaxation (characterized by Ti) and de- 
coherence (characterized by T2) of the qubits. For example, 
in the single-qubit state tomography, non-negligible T\ de- 
creases all three probabilities of the measurements, however 



non-negligible T2 reduces the probabilities of the measure- 
ments with rotations about x and y directions |28]. So the 
environmental effect on the reconstructed states is required to 
be considered in practice for more specific model. Further, 
the required quantum operations, especially two-qubit non- 
local operation, are difficult to accurately implement during 
the process of experiments. For example, the probabilities of 
theoretical calculations with qubit operations for the first and 
second qubit measurements (|l)(l|)i and (|1)(1|)2 are related 
to parameters rjj of the equivalent measurements shown in 
tables|I|and ||T| however, the measuring results of inaccurately 
experimental two-qubit operations will actually relate to not 
only these results shown in tables H] and |llj but also other 
extra terms. If these extra terms are not negligible, the re- 
constructed states might violate the properties of the positive 
semi-definiteness of the physical state p. A third error source 
is the imperfect readout of the charge qubit (for experiments, 
e.g., using single charge qubit |43], the fidelity of the read- 
out can reach 99%). The limited statistical data also affect the 
reconstruction of the states. All these imperfection can make 
the reconstructed states violate the important basic properties 
of the physical states: normalization, Hermiticity, and positiv- 
ity. In order to reconstruct a physical qubit state, in principle 
the maximum likelihood estimation of density matrices can 
be employed to minimize experimental errors. This method 
can be applied to numerically optimize the experimental data, 
which has been used in the optical systems 1 46] and more de- 
tailed discussions on this method can be referred to a very 
goodRef. JH. 

When the tomography is processed, the external flux ap- 
plied to the SQUID needs to be very quickly changed. For in- 
stance, the duration for changing $o/2 within a SQUID loop 
should at least be less than the decoherence time. Thus a pulse 
field magnetometer with a rapid sweep rate may be required 
in this experiment. If the sweep rate |48] of the pulse field 
magnetometer reaches, e.g. 10 s Oe/s, then the time to change 
<&o/2 in the loop needs about 0.25 ns for a SQUID area of 
400(/mi) 2 . 

We also notice that the number of rotations for the mea- 
sured density matrix elements to a preferable direction (e.g. y 
instead of z) grows exponentially with the number of qubits. 
How to solve this problem is still an open question. 
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